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Cold molecular ions are promising candidates in various fields ranging from precision spectroscopy 
and test of fundamental physics to ultra-cold chemistry. Control of internal and external degrees of 
freedom is a prerequisite for many of these applications. Motional ground state cooling represents 
the starting point for quantum logic-assisted internal state preparation, detection, and spectroscopy 
protocols. Robust and fast cooling is crucial to maximize the fraction of time available for the actual 
experiment. We optimize the cooling rate of ground state cooling schemes for single ^®Mg^ ions and 
sympathetic ground state cooling of ^"^MgH"^. In particular, we show that robust cooling is achieved 
by combining pulsed Raman sideband cooling with continuous quench cooling. Furthermore, we 
experimentally demonstrate an efficient strategy for ground state cooling outside the Lamb-Dicke 
regime. 


I. INTRODUCTION 

Coherent manipulation of quantum systems has been 
a goal for a wide range of fields in physics and chemistry. 
The ability to control internal and external degrees of 
freedom of atoms has enabled applications such as quan¬ 
tum computation [1-4], quantum simulation [5, 6], and 
quantum metrology [7-12]. The control over the inter¬ 
nal degrees of freedom of the atomic species is often ac¬ 
complished through coherent manipulation using laser, 
microwave or radiofrequency pulses. The external de¬ 
grees of freedom are typically controlled with various 
laser cooling techniques. Extending these coherent ma¬ 
nipulation techniques to molecular ions promises novel 
applications including multi-qubit quantum memories, 
quantum information processing [13, 14], and ultracold 
chemistry [15-17]. Moreover, molecular ions are particu¬ 
larly promising candidates for several spectroscopic tests 
of fundamental physics, ranging from a possible varia¬ 
tion of the electron-to-proton mass ratio [18-22] over an 
enhanced sensitivity to a possible dipole moment of the 
electron [23-27] to parity violation observable through 
a small energy difference in enantiomers of polyatomic 
molecules [28-31]. 

However, controlling the quantum state of molecular 
ions represents a challenge arising from their rich internal 
level structure, including rotation and vibration. As a re¬ 
sult of the absence of selection rules for vibrational tran¬ 
sitions (except for a few special cases with near-diagonal 
Franck-Condon matrices [32-34]) no closed cycling tran¬ 
sitions are available to implement laser cooling. This 
limitation can be overcome by sympathetic cooling with 
an atomic species trapped together with the molecule. 
While sympathetic Doppler cooling of molecular ions has 
been demonstrated by several groups in the past [35-39] , 
only recently sympathetic ground state cooling of CaH+ 
using Ca’*' has been achieved [79] . The approach is iden¬ 
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tical to previous successful implementations using atomic 
ions trapped in linear radio-frequency Paul traps [7, 40- 
44]. 

Besides the lowest achievable mean population of mo¬ 
tional states, h, which can reach n ~ 0.001 for a single 
ion [45], the time spent on cooling is an important aspect. 
Short cooling times reduce overhead in the experimental 
cycle and are relevant for quantum algorithms [44, 46] , as 
well as spectroscopy experiments requiring ground state 
cooling, such as quantum logic spectroscopy [7, 9] and 
photon recoil spectroscopy [11]. 

Cooling to the motional ground state of a single mode 
can be achieved via resolved sideband cooling (SBC), 
which requires the motional sidebands in the excitation 
spectrum of the system to be spectrally resolved. It has 
been demonstrated in various systems ranging from sin¬ 
gle trapped ions in Paul traps [47-50] to neutral atoms 
in optical lattices [51, 52] or optical tweezers [53] and 
micromechanical oscillators [54-56]. In the following we 
will focus on ground state cooling of trapped and local¬ 
ized one- and two-ion crystals, represented by 2-level sys¬ 
tems with two metastable electronic states (|t), [i))- The 
crystals are confined in a harmonic trap with motional 
Fock states jn), separated by the trap frequency ujt, as 
shown in Fig. 1(a). Kinetic energy is removed from the 
system by selectively driving red-sideband (RSB) transi¬ 
tions involving a change of the electronic state while re¬ 
ducing the motional quantum number, followed by spon¬ 
taneous emission back to the electronic ground state. In 
the Lamb-Dicke regime, changes in the motional state 
upon spontaneous emission are strongly suppressed [57- 
59]. 

Sideband cooling is typically implemented in two 
flavours. In continuous sideband cooling, a quench laser 
coupling the metastable excited state (j)')) to a short¬ 
lived state is applied simultaneously with the RSB laser, 
effectively broadening the linewidth of the excited state 
to optimize the cooling rate [45, 47, 60]. In pulsed side¬ 
band cooling, the RSB and the quench (also called re¬ 
pump) lasers are applied sequentially [48, 61]. 

The RSB coupling strength depends strongly on the 
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motional state [57, 62] and outside the Lamb-Dicke 
regime even exhibits points of vanishing coupling for cer¬ 
tain initial motional states, effectively disabling further 
cooling beyond these points (see Fig. 2). While this 
regime has been studied theoretically [63-65], no experi¬ 
mental investigations are known to us. 

Here, we demonstrate fast and robust sympathetic 
ground state cooling of a molecular ion along 

one direction of motion using a ^®Mg“'' cooling ion. Start¬ 
ing with a Doppler-cooled single trapped ^^Mg"'' ion, we 
investigate in Sec. II a novel repumping scheme for pulsed 
Raman sideband cooling, in which the excited electronic 
state in the cooling cycle is effectively quenched to the 
electronic ground state via coupling to a short-lived ex¬ 
cited state. We demonstrate that the requirements on 
meeting the optimum RSB pulse length are significantly 
relaxed for the new quasi-continuous scheme compared 
to conventional pulsed Raman SBC. As a consequence 
of the large linewidth of the Doppler cooling transition 
in ^^Mg’*', a significant amount of motional state popu¬ 
lation is trapped above the point of vanishing coupling 
strength for RSBs. We employ second order RSB tran¬ 
sitions to sweep the population beyond this point and 
demonstrate in Sec. IV that ground state cooling with 
Lamb-Dicke factors as large as 0.45 and motional level 
populations up to n ~ 120 becomes possible by employ¬ 
ing RSBs up to 8 **' order. By optimizing the ground state 
cooling rate in terms of pulse lengths for 1 ®* and 2 "'^ order 
RSB and the time spent on cooling with each sideband 
order, a total cooling time as short as 500 ps for cooling 
a single ion from Doppler temperature (h ~ 10) to the 
motional ground state is demonstrated. The experimen¬ 
tal results are supported by numerical Master equation 
simulations of the system (Appendix A). In Sec. Ill this 
scheme is extended to cool a ^^MgH’^/^^Mg^ ion crystal 
to the ground state by interleaved cooling of both axial 
modes. After optimizing the RSB pulse length for each 
mode and the time spent on cooling each, a total cooling 
time of 2.5 ms is achieved, resulting in a mean residual 
motional excitation of 0.06(3) and 0.03(3) for the in- and 
out-of-phase mode, respectively. 


II. SBC OF A SINGLE ION 
A. Experimental setup 

The current work is based on previous results [ 66 ], 
where we demonstrated a pulsed Raman sideband cooling 
scheme using a single laser system for ^^Mg"*". In brief, 
a 25Mg+ ion is trapped via isotope-selective photoioniza¬ 
tion in a linear Paul trap with axial and radial motional 
frequencies of 2.21 MHz and ~ 5 MHz, respectively. The 
relevant levels for SBC of a single ^^Mg“'' ion are shown in 
Fig. 1(b). A frequency-quadrupled fiber laser ("P 3/2 laser 
in the following) provides the radiation for Doppler cool¬ 
ing, coherent control, and state detection on the ^®Mg“'' 
ion. Doppler cooling on the [^Si/ 2 , F = 3, ttt-f = 3) = [),) 
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Figure 1. (Color online) Principle of resolved sideband 
cooling and implementation in (a) Resolved 

SBC requires the linewidth F of the transition to be smaller 
than the trap frequency ljt, so that individual transitions 
can be selectively driven. With the cooling laser tuned to the 
1*** order RSB, one quantum of motion is removed from the 
system in each absorption event. In the Lamb-Dicke regime, 
a change in motional state is suppressed upon decay back to 
the [4) state, (b) Relevant level structure for SBC of a 
ion. Doppler cooling and Raman transitions are performed by 
coupling the ^Si /2 and ^P 3/2 states. Compared to a previous 
implementation of SBC [66], a second laser is added for re¬ 
pumping the [aux) state via the ^Si /2 ^ ^Pi /2 transition. RF 
radiation couples the states Id) and |aux). 


O pP 3 / 2 ,F’ = 4:,mF = 4 ) transition for 400ps yields a 
fi ^ 10 as the starting point for the ground state cooling 
sequence. Coherent electronic and motional state ma¬ 
nipulation is implemented via Raman transitions, cou¬ 
pling the 14) and |t) = 1 ^ 81 / 2 ,!^ = 2,mF = 2) states. No 
change in the motional state correspond to carrier (CAR) 
transitions, while addition of phonons correspond to blue 
sideband (BSB) transitions upon changing the electronic 
state from |4) to [t)- The Raman beams are generated 
via acousto-optic modulators from the ^P 3/2 laser and 
are detuned with respect to the ^P 3/2 state by around 
9.2 GHz. The beams have tt and polarization to maxi¬ 
mize the coupling strength and are at right angle to each 
other with the difference wavevector aligned along the 
axial direction of the trap. This results in a Lamb-Dicke 
parameter of 77 = 0.3 [66] for the axial mode of interest 
here. The limited detuning of the Raman laser beams 
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Figure 2. (Color online) Motional state population and 
effective Rabi frequencies, (a) Typical population distri¬ 
bution over motional states after Doppler cooling, correspond¬ 
ing to n « 10. (b) Effective Rabi frequencies as a function of 
the motional state calculated with a Lamb-Dicke parameter 
77 = 0.30 according to Eq. (1). The effective Rabi frequencies 
for different RSB orders have zero points at different motional 
excitations, solid/dashed/dotted-. order RSB. 


leads to off-resonant excitations, which result in dephas¬ 
ing and population loss during coherent manipulation as 
further discussed in Sec. IID. This effect limits the detec¬ 
tion fidelity of the motional ground state population and 
the final n detectable after cooling. It can be reduced 
by employing a separate Raman laser system with larger 
detuning. 

In the previously implemented SBC scheme [66] the 
Rabi frequency ^n',n of the Raman transition between 
motional state jn) and |n') was calculated according to 
[57, 62] 

= floexp (1) 

where Oq is the CAR Rabi frequency, (n^) denotes 
the smaller (larger) of n, n', and T“(x) are the gener¬ 
alized Laguerre polynomials. This allowed the deter¬ 
mination of the TT-time for complete population trans¬ 
fer on the 1®* and 2"'^ order RSBs (see Fig. 2). This 
way, RSB 7r-pulses were applied to sweep the popula¬ 
tion starting from n = 40 to the motional ground state. 
After each RSB pulse, multiple repump cycles involv¬ 
ing optical excitation of the |t) —>■ ^P 3/2 transition and 
RF state inversion pulses between ]aux) and |t) are nec¬ 
essary to transfer the population back to the ]),) state, 
since the excited state in ^P 3/2 can also decay into the 
]aux) = \^Si/ 2 ,F = 3,mi? = 2) state. The cooling speed 
of this implementation was mainly limited by the radio 
frequency (RF) coupling strength and the required num¬ 
ber of sweeps (2-3) to reach the ground state. 


B. Cooling scheme 

In the approach described in the following, we signifi¬ 
cantly reduce the cooling time by adding a dedicated re¬ 
pump laser for faster repumping of the jaux) state via the 


I^Pi/ 2 ,P" = 3,mi? = 3 ) state and by using RSB pulses 
with equal length. The experimental sequence starts with 
Doppler cooling. For SBC, we first apply a series of 2“'^ 
order RSB pulses followed by a series of 1®* order RSB 
pulses of fixed lengths. After each SBC pulse a short 
(tr = 3 fj.s) optical repumping pulse using the beam 
of the Raman laser tuned to resonance and the ^Pi /2 
laser is applied to clear out the population left in jaux) 
and ]'(). In the experiment the actual elapsed time be¬ 
tween two SBC pulses equals = 5 jis caused by delays 
from the control electronics. In addition, we can use the 
^Pi /2 laser and the RF coupling between jaux) and jf) 
to implement a quench coupling during the Raman RSB 
pulses. This quench coupling opens up an additional de¬ 
cay channel for the |t) state back to the ground state 
I),) with a controllable decay rate, implementing a fusion 
between pulsed and continuous sideband cooling. After 
SBC, the population in the motionally excited states is 
probed by driving a RF 7r-pulse to transfer all population 
from I),) to It), followed by a stimulated Raman adiabatic 
passage (STIRAP) pulse on the 1®* blue sideband (BSB) 
]67]. This maps motionally excited states (n > 0) onto 
the state ]t), while the motional ground state population 
(n = 0) remains in the state |t)]68]. 

We optimize the cooling rate for a fixed total cooling 
time Tc by changing the amount of time spent on cool¬ 
ing via 2"'^ and 1®* order RSBs, characterized by a time 
scaling factor a. Depending on the pulse length of the 
pt ( 2 'Ki) order RSB pulses tm {tm), we apply Nb .2 2"'^ 
order and iVm 1®* order RSB pulses 


A^R2 = 


A'ri = 


aTc 
tR2 _ 

(1 - a)T, 


( 2 ) 


where the sign [xj = max{n € Z, n < a:} denotes the 
floor function. Short padding pulses are added to keep 
fixed. The total time spent on ground state cooling dur¬ 
ing the cooling sequence as shown in Fig. 3 is expressed 
as 


Ttotal — Tc + Nc ■ t[., (3) 

where W is the total number of repump pulses. 

Since we are interested in the cooling time, we inves¬ 
tigate the dynamics of SBC by probing the population 
in the motional ground state with the STIRAP sideband 
pulse for different SBC times Tc instead of optimizing for 
lowest h. Assuming a constant cooling rate W during the 
SBC cycle, the mean occupation of the motional states 
decays exponentially [69] as 

h(t) = hf(l — e"’^*)-|-(4) 

where hi and hf are the initial and final mean occupa¬ 
tion, respectively, and W is the cooling rate. Assuming 
a thermal distribution with mean occupation of n(t) af¬ 
ter SBC for a duration t, the population in the motional 
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Figure 3. (Color online) Sequence for single ion sideband 
cooling. Sequence for SBC a single The sequence 

starts by repeating Nia. 2“"^ order RSBs, followed by A^ri 1®* 
order RSBs with fixed pulse lengths. After each SBC pulses, 
repumping pulses consisting of ^Pi /2 and a beams are applied 
to bring the ion back to |4,) state. 
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ground state can be expressed as 
^ 1 + ntil - 

I (5) 

1 + nf + (fii — n{)e~*/'^o 

The desired cooling time constants Tq = 1/W for sets 
of parameters , 1^,2 > <x) are extracted by fitting data 
to this model function in which only one parameter is 
changed and hi and hf are common fit parameters. By 
choosing Tc ^ 7 ■ Tq, we ensure reaching the steady state. 


C. Optimization of experimental parameters 

Using the procedure described in the previous section, 
we vary the pulse length of the RSB pulses, the time scal¬ 
ing factor a and the optical power of the ^Pi /2 repump 
laser to optimize the cooling rate. Furthermore, we ex¬ 
plore the transition between pulsed and quench cooling 
by adding the ^Pi /2 repump laser and RF coupling be¬ 
tween the |aux) and |t) states during the RSB pulses. 
A numerical simulation based on optical Bloch equations 
described in Appendix A supports our findings quantita¬ 
tively. 


1. Pulse length optimization 

We first optimize the cooling time constant as a func¬ 
tion of the pulse length for the 1®* and 2"'^ order RSB 
pulses. For each parameter set, a scan of the SBC time Tc 
is performed, which allows us to derive the corresponding 
cooling time constant as shown in Fig. 4(d). In Fig. 4(a) 
and 4(b), the experimentally determined cooling time 
constants Tq are plotted against the pulse lengths. The 
shortest Tq is obtained for 1®* and 2"'^ order RSB pulse 


Figure 4. (Color online) Dependence of the cooling time 
constants on experimental parameters. Cooling time 
constants as a function of the pulse lengths of the (a) 2 "'^ order 
and (b) order RSB pulses with a = 0.5 while fixing the 
non-scanned pulse length near its optimal value, blue point 
(red square)-. Experimentally determined cooling constants 
with the quench coupling on (off) during the RSB pulse. The 
quench coupling induces an effective decay rate of 42 ms^^ 
from the |4) state to the If) state, (c): Cooling time constant 
as a function of the time scaling factor with tB .2 = 10 ps and 
till = 10 ps. The lines in (a)-(c) are the result of Master 
equation simulations using the experimental parameters, (d) 
Residual motional excitation as a function of the total cooling 
time Tc. The line is the fit to the experimental data, which 
gives the cooling constant Tq. 


lengths of around 10 ps. The experimentally determined 
TT-time for a RSB pulse starting at \n = 1) is « 16 ps. 
Averaging over the 7r-times for 1®* and 2'^'^ order RSBs ac¬ 
cording to Eq. (1), weighted by the thermal occupation of 
the initial states, limited for each sideband to the states 
where the coupling strength of the respective sideband 
dominates, gives for both sidebands average 7r-times of 
^10 ps in agreement with the experimentally optimized 
values. The fact that a constant 7r-time is sufficient for 
efficient cooling is a consequence of the relatively small 
variation of the Rabi frequencies in the relevant range of 
motional levels as can be seen from Fig. 2. 

In standard pulsed Raman SBC, the cooling time con¬ 
stant strongly rises upon deviation from the optimum 
pulse lengths [red points and dashed lines in Fig. 4(a) 
and (b)]. This dependence becomes much weaker for long 
pulses [blue points and solid lines in Fig. 4(a) and (b)[ 
by applying the quench coupling during the RSB pulses 
which opens a decay channel for the [t) state. In the ex¬ 
periment, this decay rate was adjusted to about 42 ms^^ 
or a decay constant of about 24 ps to achieve the highest 
cooling rate. At this decay rate, the effect of repumping 
is negligible for pulses shorter than the optimum pulse 
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Figure 5. (Color online) Rabi excitation of after 

sideband cooling, (a-c) Frequency scans of Raman transi¬ 
tions over RSB (red), carrier (black) and BSB (blue) 
transitions after SBC of a single 


length and the scheme is equivalent to standard pulsed 
SBC. For RSB pulses longer than the optimum pulse 
length, the ion in the |t) state decays back to the jt) 
state through the new channel. The ion cycles between 
the 14) and the 14) states and therefore more than one 
phonon can be removed within a single RSB pulse. The 
scheme thus becomes equivalent to continuous SBC. This 
quasi-continuous cooling scheme is insensitive to the ex¬ 
act pulse length of the RSB pulses and provides high 
robustness against intensity/pointing fluctuation of the 
Raman lasers. 


2. Time scaling factor optimization 

The time scaling factor a as defined in Eq. (2) dis¬ 
tributes the SBC time Tc into the time spent on the 
order RSB Trz = aT^ and on the order RSB 
Tm = (1 — a)Tc. Due to the dependence of the Rabi fre¬ 
quency on the motional quantum number (see Fig. 2 and 
Eq. (1)), the 2"'^ order RSB pulses are more efficient in 
cooling the population in motional states n > 20 at the 
starting stage of the SBC cycle. For the lower motional 
states, the 1®* order RSB pulses become more efficient, so 
that an a that is too large also increases the SBC cooling 
time constant (Fig. 4c). Since the variation of 1®* and 
2 nd RSB Rabi frequencies with n are very slow, 

the cooling time spent on either order only significantly 
influences the cooling time constant for the extreme cases 
of a ^ 0 and a —>■ 1 and remains otherwise flat. 


D. Cooling results 

The final fi after SBC is extracted from the red and 
blue sideband excitation. Assuming that the ion stays in 
a thermal distribution after SBC, the ratio of excitations 
on 1®* red and blue sideband fulfills the following relation 
[70] 


p''®‘’(t) n 
~ l + n 


( 6 ) 


where p^^'°(t) and p^^'°(t) are the excitation probabilities 
on the 1®* RSB and BSB at time t, respectively. For this 
analysis we subtracted a constant signal background from 
the data. With the SBC sequence described in Sec. IIB, 
we achieve n = p'^^'°/(p'°^'° — p^^'°) ~ 0.01(2) as shown in 
Fig. 5(a)-(c) after a cooling time of Tc = 500 ps. This 
represents a reduction in cooling time by more than a 
factor of 30 compared to reference [66]. 


III. SYMPATHETIC GROUND-STATE 
COOLING OF A MOLEGULAR ION 


In this Section we adapt the quasi-continuous cooling 
scheme for sympathetic cooling of a molecular ion using 
25Mg+ as the cooling ion species. 


A. Loading and sympathetic Doppler cooling of 
^^MgH+ 

A two-ion crystal consisting of a ^®Mg“'' and a ^^MgH"'' 
ion is prepared by isotope-selective photo-ionization load¬ 
ing [71] of a pair of ^^Mg+ and ions with the 

Doppler cooling laser tuned near the resonance of the 
25Mg+ cooling transition. Then, the laser is tuned near 
to the cooling resonance of ^'^Mg'*' and hydrogen gas is 
leaked into the vacuum system through a leak valve, in¬ 
creasing the pressure up to ^ 5 x 10“® mbar. After a 
photo-chemical reaction during which the excited ^^Mg“'' 
ion reacts with a hydrogen molecule to form ^"'^MgH''' 
[72], the fluorescence of the ^^Mg“'' ion vanishes. After 
closing the leak valve and tuning the laser back to the 
25Mg+ resonance, cooling commences and the mass of 
the dark ion is determined via mass spectrometry using 
parametric heating of the two-ion crystal [73]. Since the 
mass of the two ions is almost identical, we do not expect 
significant deviations from Doppler cooling temperature, 
even in the presence of additional heating [74]. 


B. Two-mode sympathetic ground-state cooling 
sequence 

For sympathetic ground state cooling of the molecu¬ 
lar ion we use a slightly modified pulse sequence com¬ 
pared to Sec. IIB. In contrast to the case of a sin¬ 
gle ^®Mg“'', the motion of the ions is described by two 
modes along the axial direction, the in-phase (ip) and 
the out-of-phase (op) mode with secular frequencies of 
= 27r X 2.21MHz and = 27r x 3.85 MHz. The 
Lamb-Dicke parameters for the coupling of the Raman 
lasers to the ^^Mg^ ion are pip = 0.21 and pop = 0.16. 
With these Lamb-Dicke parameters, the effective Rabi 
frequencies show no zero-crossings over the range of trap 
levels with significant population after Doppler cooling, 
so no higher order RSB pulses are necessary. However, 
with Lamb-Dicke factors as large as these, cooling of a 
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Figure 6. (Color online) Sequence for sympathetic side¬ 
band cooling of a two-ion crystal. The in-phase and 
out-of-phase axial modes of a 2-ion crystal 

are cooled in an interleaved fashion to ensure simultaneous 
cooling. 


single mode is not sufficient to enable high fidelity op¬ 
erations involving only this mode. The other mode acts 
as a spectator mode which modifies the effective Rabi 
frequency of the mode of interest depending on its mo¬ 
tional state [62]. Therefore, we have implemented an 
interleaved pulse sequence for SBC both axial modes si¬ 
multaneously as shown in Fig. 6. Depending on the time 
scaling factor a', which in this case distributes the total 
SBC time into time spent for RSB on the ip mode 
(1 — a')Tc and time for RSB on the op mode a'T^, we 
apply 


Nip = 
Nop = 


(1 - a')T, 

tip 



(7) 

( 8 ) 


pulses on the ip and the op mode, respectively. Here 
tip, top are the pulse lengths of the RSB pulses on the ip 
and the op mode. In the case Nip > Nop {Nip < Nop) 
we start the SBC cycle with iV^es = |tVip — A^^pj pulses 
on the ip (op) mode, followed by Nop {Nip) pulses on 
the ip and the op mode in an interleaved fashion. After 
every single SBC pulse, a repumping pulse as described 
above is applied to clear out the |t) and jaux) states in 
addition to the quench coupling present also during the 
RSB pulses. Optimizations similar to the case of a single 
25Mg+ are performed to minimize the total duration of 
SBC. 


C. Pulse length optimization 

The length of the RSB pulse on the ip (op) mode is 
optimized by fixing the other mode’s RSB pulse length. 
As described in Sec. IIB, the cooling time constant is 
derived from the motional ground state population mea¬ 
sured by performing a STIRAP pulse resonant with the 
P* BSB of either mode for different cooling times To- We 
spent an equal amount of time for cooling on both modes 



Figure 7. (Color online) Dependence of the cooling time 
constant on experimental parameters for sympathetic 
sideband cooling a two-ion crystal, (a) The cooling time 
constant as a function of the optimal pulse length of the Ra¬ 
man RSB pulses on the ip {blue circles) and op {red squares) 
mode, (b) The cooling constant as a function of the optimal 
time scaling factor a'. The lines connect the points and are 
guides to the eye. 


{a' = 0.5) and ensure near steady-state conditions by 
performing cooling up to Tc ~ 2.5 ms. The highest cool¬ 
ing rate is observed with a pulse length of 15 /is (20 ps) 
for the RSB on the ip (op) mode as shown in Fig. 7a. 
With a TT-time of « 21 p,s and k, 26.5 ps for 

the 1®* RSB starting in the |n = 1) state, the Rabi fre¬ 
quency according to Eq. (1) averaged over the thermal 
occupation of motional states results in a mean 7r-time 
of 14 ps and 25 ps for the ip and op mode, respectively. 

As in the single-ion case discussed in Sec. IIC, the 
increase in the cooling time constant for pulse lengths 
longer than the optimum pulse length is rather mild as a 
consequence of the quench coupling during RSB pulses. 


D. Time scaling factor optimization 

The time scaling factor in the SBC sequence for cool¬ 
ing a two-ion crystal divides the total cooling time To into 
the time spent on cooling the ip mode and the op mode 
for the optimized pulse lengths derived from the mea¬ 
surements in the previous section. The measured SBC 
cooling time constants for the two modes as a function of 
a' are shown in Fig. 7b. Although an a' that is too small 
(large) accelerates the cooling of the ip (op) mode, the 
cooling time of the opposite mode increases. The optimal 
value for the time scaling factor a' ps 0.5 is determined 
by the point where SBC of both the ip and op mode are 
achieved with the shortest total duration of the SBC cy¬ 
cle. Interestingly, this corresponds to a different number 
of RSB pulses for the two modes. 


E. Cooling results 

Fig. 8 shows frequency scans of the carrier and both 
RSB and BSB transitions after SBC the ip and op modes 
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Figure 8. (Color online) Rabi excitation of after 2-mode SBC a crystal. Frequency scans of 

Raman transitions over the RSB (red) of the (a) op and (b) ip modes, (c) the carrier (black), and the F** BSB (blue) (d) 
op and (d) op modes after SBC both modes on ^®Mg+. 


along the axial direction with optimized parameters. The 
final h is determined from the red and blue sideband 
excitations of each mode as described in Sec. IID. As 
before, the offset from off-resonant Raman excitations 
has been subtracted for the analysis. We reach a hip = 
0.06(3) for the ip mode and hop = 0.03(3) for the op mode 
after a total cooling time Ttotai ~ 2.5 ms. Compared to 
an extension of the SBC scheme described in reference 
[ 66 ] to 2 -mode ground state cooling, the optimized new 
scheme reduces the cooling time by a factor of 8 . 


IV. SIDEBAND COOLING BEYOND THE 
LAMB-DICKE REGIME 

For some applications large Lamb-Dicke factors are de¬ 
sirable, since they enhance the sensitivity of the ion’s 
motion to small forces. This is the case for photon 
recoil spectroscopy [ 11 ] or the non-destructive internal 
state detection of a molecular ion using oscillating dipole 
forces [75-77]. In the following we extend the single¬ 
ion ground state cooling scheme presented in Sec. II to 
enable efficient ground state cooling outside the Lamb- 
Dicke regime, which is readily extended to the multi-ion 
case discussed Sec. III. 

At a trap frequency of 1 MHz and a temperature of 
1 mK theoretically achievable with Doppler cooling a 
25Mg+ ion, states up to n ^ 120 are significantly popu¬ 
lated, leaving less than 0.3% population in levels n > 120 
(Fig. 9a). At this trap frequency, the Lamb-Dicke param¬ 
eter in our system is 77 = 0.45. The effective Rabi fre¬ 
quency for RSB transitions depends strongly on the trap 
levels and shows several points of vanishing coupling over 
the range of trap levels with significant population (see 
Fig. 9(b) and Eq. (1)). Employing pulses on the 2"'^ or¬ 
der sideband as in the scheme for a single ^®Mg“'' is no 
longer sufficient. Instead, we employ as many higher or¬ 
der sidebands as necessary. For larger motional states, 
successively higher order sidebands exhibit a maximum 
in their coupling rate. Moreover, pulses on higher order 
sidebands are more efficient since more than one phonon 
is removed per sideband pulse. For simplicity, we split 
the total cooling time Tc equally between all the side¬ 


band orders and apply to all of them the same pulse 
length, thus extending the sequence shown in Fig. 3 to 
higher-order modes. This is a valid approach, since the 
variation in maximum Rabi frequency across the higher 
order sidebands is small. Simultaneously to each RSB 
pulse, we apply the quench coupling to further reduce 
the sensitivity to the optimum pulse length. After each 
RSB pulse, a repumping pulse as described before is ap¬ 
plied. The optimal number of applied sideband orders, 
derived from measurements of the cooling time constant 
as a function of the maximally applied sideband order, is 
shown in Fig. 10a. In the optimal sequence, we need to 
apply sideband pulses up to the 8 *^ order, which confirms 
the prediction of Fig. 9b. 

More generally, with known Lamb-Dicke parame¬ 
ter and thermal distribution over motional states after 
Doppler cooling, a corresponding SBC strategy can be 
used to achieve the most efficient cooling as illustrated 
in Fig. 10b, where we define the cooling efficiency as 
/3(n) • For a selected upper trap level to be ad¬ 

dressed by cooling and a known Lamb-Dicke parameter. 



Trap levels 


Trap levels 


Figure 9. (Color online) Motional state population and 
effective Rabi frequency for small trap frequency, (a) 

Distribution of the motional levels at a motional frequency of 
ujt = 1 MHz and a temperature of T = 1 mK corresponding 
to the Doppler cooling limit of the ^^Mg’*' ion. Less than 
0.3% of the population are in the trap levels \n > 120). (b) 
Effective Rabi frequencies at 77 = 0.45 for different RSB orders 
are shown as a function of the trap levels. For ions located in 
higher trap levels, RSB pulses of higher orders /3(n) are more 
efficient for cooling. 
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Figure 10. (Color online) Determination of the maxi¬ 
mum sideband order, (a) The cooling time constant as a 
function of applied sideband orders. The optimal number for 
the sideband orders determined experimentally confirms the 
prediction shown in Fig. 9(b). (b) RSB orders with highest 
effective cooling rate as a function of the trap level and the 
Lamb-Dicke parameter. With known Lamb-Dicke parameter 
and motional distribution a corresponding SBC sequence can 
be adopted. 


forces through the excitation of the ion’s motion from 
the ground state, such as in photon recoil spectroscopy 
[11] and the detection of electric fields [61, 78[. 

The demonstrated motional ground state cooling of a 
molecular ion, in particular when combined with cooling 
at low trap frequency, represents an important step to¬ 
wards the implementation of non-destructive state prepa¬ 
ration and detection techniques [75-77] based on the de¬ 
tection of small optically-induced and state-selective dis¬ 
placement forces acting on a trapped molecular ion. 
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the sideband order with highest cooling efficiency is plot¬ 
ted. This strategy can be combined with multi-mode 
cooling through interleaved pulse sequences as shown in 
Sec. III. 


V. SUMMARY & DISCUSSION 

A quasi-continuous sideband cooling scheme was in¬ 
troduced by applying a quench coupling to the excited 
state in pulsed Raman sideband cooling. For long Raman 
pulses this scheme allows multiple RSB-spontaneous 
emission cycles, reminiscent of continuous quench cool¬ 
ing. This feature makes SBC more robust and signif¬ 
icantly reduces the optimization required for optimum 
cooling, since RSB pulses that are longer than the 7r-time 
for a specific initial motional state contribute more effi¬ 
ciently to cooling. Using this scheme, we performed an 
optimization of the cooling rate of single-ion ground state 
cooling and two-ion sympathetic ground state cooling, 
significantly outperforming previous results [66[. These 
findings are important to minimize the time required for 
ground state cooling, thus reducing overhead in experi¬ 
ments with trapped ions. The results are applicable to 
other commonly used hyperfine qubit systems such as 
^Be+j ^^^Cd+, or ^^^Yb+. A variation of the scheme can 
be applied to other systems like ^°Ca+, in high 

magnetic field or optically trapped neutral atoms. 

We have also demonstrated experimentally that SBC 
to the ground state is possible outside the Lamb-Dicke 
regime by employing higher-order RSB transitions, thus 
confirming theoretical predictions [63-65] . This regime 
is particularly relevant for systems where the Doppler 
cooling linewidth significantly exceeds the motional fre¬ 
quency as is the case for optically trapped neutral atoms, 
or in situations where the Lamb-Dicke factor is delib¬ 
erately chosen as large as possible. This is the case 
for example in experiments which aim to detect small 


Appendix A: Numerical simulation for SBC 

The dynamics of the system during quasi-continuous 
SBC is modelled using optical Bloch equations, where 
two electronic states |t) and [J,) and 80 trap levels 
(Fig. 11) are considered. The Raman RSB pulses 
are included as a resonant coupling between j),, n) and 
It, n — P) with P indicating the sideband order used. The 
spontaneous decay during the Raman cooling cycle via 
the auxiliary states [80] is implemented with an effec¬ 
tive decay rate Fefi determined experimentally. This ap¬ 
proach neglects effects from the actual multi-level elec¬ 
tronic structure of ^^Mg'*'. Repumping after applica¬ 
tion of the RSB is implemented through a projection 
of the population onto the electronic ground state jj,). 



n=0) |n=1> |n=2) |n=3) |n=4) 


Figure 11. (Color online) Relevant levels for numerical 
simulation of the dynamics during the Raman SBC 
cycles. The Raman lasers are tuned to be resonant with 
the order RSB (/3 = 1 in graph). Through the dissi¬ 
pative channels either from off-resonant excitation or via a 
repumping process followed by spontaneous emission, the ion 
is reinitialized to the [f,) with a decay rate of 0.005 s^^ and 
0.47 s^'^, respectively. Spontaneous emission happens on the 
carrier transition with a probability 1 — ^ and the remaining 
small fraction ^ happens on the RSB transitions. 
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The spatially averaged Lamb-Dicke factor for sponta¬ 
neous emission along the axial direction is fj = 0.134. 
Heating from spontaneous emission is taken into account 
through a branching ratio of 1 — ^ ^ for emission on 

the CAR and RSB transition. Here ^ is on the order 
of 3 • 7^ ~ 0.05, where the factor 3 considers the multi¬ 
ple scattering events until the ion falls into the |4.) state. 
Combining all the ingredients above, we end up with the 
following optical Bloch equations 

P4,n,4,n = i ^ [P4,n,t(n-/3) ~ Pt(n-;8),4.n] 

-I- (1 — • Tefiptn.tn 

+ ? • reffPt(n-l),t(n-l) (Al) 

Pin,t(n-/3) = ~ P^{n-P),in 

.^n—/3,n r l 

— * 2 “ Pt(n-;S),t(n-,fl)J 

-^P4.n.t(n-/3) (A2) 

Pt(n-/3),t(n-/3) = * ” 2 ^’” [^t(n-/3).tn “ Pgn,t(n-/3)] 

- reflPt(n-/3),t(n-/3) (A3) 

with the density matrix elements defined in the usual 
way, e.g. 

Pin,t(n-/3) = (i, n] p It, (n - /?)) (A4) 


and as the effective Rabi frequency coupling the 

two states |t, n) and |t, n — 13). This system of differential 
equations is solved numerically to produce the evolution 
during a single Raman cooling pulse. 

For every set of parameters (Tc, a, tm and tR 2 ), the 
evolution during the pulse sequence as shown in Fig. 3 
is computed by repeatedly solving the equations above. 
The actual evolution time for the atomic system dur¬ 
ing each pulse is reduced by 1 ps from the chosen pulse 
lengths to include realistic experimental pulse area re¬ 
duction. The final population in the motional ground 
state is considered as the signal detected by the STIRAP 
pulse and corrected by an amplitude a = 0.7 and an off¬ 
set b = 0.23 reduction due to experimental imperfections 
according to 


y — a X (ytheory ^)- (A 5 ) 


A scan over the SBC cooling time Tc reproduces the ex¬ 
perimental data, which are processed in a similar way as 
described in Sec. HB. 
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